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Specific Heat of a Three Dimensional Metal Near a T = 
Magnetic Transition with Dynamic Exponent 2; = 2, 3, or 4 
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Abstract 

We derive expressions for the universal contribution to the specific heat of 
a three-dimensional metal near a zero-temperature phase transition with dy- 



^ ■ namic exponent z = 2, 3, or 4. The results allow a quantitative comparison of 



theory to data. We illustrate the application of our results by analyzing data 
for Cei_xLuxCu2Si2, which has been claimed to be near a quantum critical 



G^ ■ point. 
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I. INTRODUCTION 

Many strongly-correlated metallic compounds |l]] exhibit " non-fermi-liquid" behavior in 
the sense that the low temperature (T) specific heat (C) exhibits more rapid temperature 
dependence than expected from the standard expression for a fermi liquid with fermi tem- 
perature Tp, which is 

T = 7>^ + «i?'°«^ + °(7f ) + ■■■] (") 

with A = ^^^- n (n denotes the electronic density) and B given in terms of Landau pa- 
rameters P]. Many authors have suggested that the non-fermi-liquid behavior is due to the 
proximity of the non-fermi-liquid materials to a zero temperature critical point f^-^ . One 
class of proposed critical points would separate a T = magnetically ordered phase from 
a phase with no long range order |0-^; another class involves critical points occuring in 
models of isolated impurities 0. To experimentally distinguish these intriguing suggestions 
it is necessary to compare the data to theory. In this communication we provide quantitative 
expressions for the universal contributions to the specific heat of a three dimensional metal 
near a zero temperature magnetic — non-magnetic phase transition with a dynamical expo- 
nent z = 2,3, or 4. The case z = 2 describes antiferromagnetic transitions not driven by 2fc^ 
instabilities or nesting; the case z = 3 describes ferromagnetic transitions in a clean metal 
and also describes the recently discussed "gauge theories" of strongly correlated systems; 
the case z = 4 would describe a dirty ferromagnet in the temperature regime in which the 
randomness leads to spin diffusion but does not otherwise change the critical properties. 

Our results are obtained by solving scaling equations derived previously 0. The scal- 
ing equations imply that the critical behavior is in the Gaussian universality class, with a 
dangerously irrelevant operator. The results depend on three dimensionless parameters: re- 
duced temperature t = ^ (where T* is a microscopic temperature presumably of the order 
of the fermi temperature Tp), a control parameter r (which depends on temperature and 
on a Hamiltonian parameter such as pressure or doping whose variation tunes the material 



through the T = transition, which is defined to occur at r(T = 0) = 0) and u, the coef- 
ficient of the " dangerously irrelevant operator" . The results of course also depend on the 
spatial dimensionality d and the dynamical exponent z, and the results we write are only 
valid for d + z > A. A qualitative phase diagram is shown in Fig.l. The axes are reduced 
temperature , t, and the zero-temperature control parameter 6 = r{T = 0). A phase transi- 
tion line tc{S) separates a magnetically ordered phase (shaded region) from a phase with no 
long range order. We give results for the specific heat coefficient 7 = C/T in the disordered 
phase. The behavior in the ordered phase depends on details of the order. Too many special 
cases arise to discuss here. In the disordered phase we find that the singular part, •jsing, of 
7 = ^ may be written 

ls^ng=Y;t^ Qd^z \-^,, I (1-2) 

where gd,z is a universal scaling function which we calculate. The leading behavior of '')sing is 
given by |^t^ gd,z i]^^ j and may be computed directly from the Gaussian model. The 
function gd,z{-[iji,0) describes the crossover between a low-temperature (so-called quantum) 
regime with t <^ ra and a high-temperature (so-called classical) regime defined by t ^ ra. 
A generally valid expression for Qs^zix, 0) is given in Eq. (P3|). Expansions for the quantum 
and classical limits are given in Eqs. ( p.l|) -( ^^ . We also investigate the dependence of 
9d,z{x,y) on its second argument and show that it leads to "corrections to scaling", which 
are then given as a power series in u (times appropriate combinations of t and r). The 
only important correction to scaling arises from the nonanalyticity in gd^z{x,y) occurring at 
y = oo, i.e. r = 0. The equation r(t,6,u) = defines the true transition curve tc{S). A 
precise expression is given in Eq. ( p.7|) . To an extremely good approximation, this transition 
is a second order mean field transition and is characterized by a specific heat jump which 
we calculate. The result is Eq. ( |3.8| ). We also calculate the Gaussian fiuctuation correction 
to this specific heat jump. This is given in Eq. ( |3.9| ). 

The rest of this paper is organized as follows. In Section |I| we write down the model to 
be solved and indicate the method of solution. In Section |T| we give the results for d = 3 



and z = 2,3,4. In Section |^ we illustrate the comparison of the results to data [g by- 
comparing (73 2 to data on Cei_2,Lu2,Cu2Si2, a three-dimensional material which is argued to 
be near a T = antiferromagnetic transition. Section [V| is a brief conclusion. Details of 
various calculations are given in Appendices. 

II. FORMALISM 

We follow the previous approach [^|10| of integrating out the fermionic degrees of freedom 



and writing the partition function as a functional integral over an N-component bosonic 
order parameter which represents the ordering field. This procedure is not useful for 
superconducting transitions, transitions at nesting wavevectors of nested fermi surfaces, or 
indeed transitions at wavevectors Q connecting two points on the fermi surface with parallel 
tangents (e.g. Q = 2kp for a spherical fermi surface) because the resulting bosonic action 
is singular. For further details see Refs. 0,y,|T^. Recently, Sachdev |Tl[| has questioned 



whether the action used in Refs. 0,0 is adequate to describe magnetic phase transitions 
in which the order parameter has 0(3) or higher symmetry because of the neglect of terms 
giving rise to the precession terms in the spin-wave equation of motion, but it has been 
recently argued that the action of Refs. [pi[TO| is in fact correct even in this case [|l2l . We do 



not discuss this issue here, but assume the validity of the approach of Refs. [I,!!^], which is 
certainly correct for low symmetry order parameters and may apply to higher symmetries. 
The partition function is 

Z= fV(fexp{~S[$]) (2.1) 

— * 

and the action iS[0] is given by 

S[^] = 5(°) + 5(2) [0] + S^'^ [0] + . . . (2.2) 



(27r; 



So + i^okY' '^" 



T(^o^) 



Mk)(p^n{-k) (2.3) 



n '' 

Sm = uoV' E /'^|^---(|%'^^n.(fci)...0„.(^4)5(E4)5E.- (2.4) 



ni,...,n4 



Here S^^^ is the logarithm of the contribution to the partition function from the degrees of 
freedom which have been integrated out. It depends on T and other external parameters. 
The parameter 60 tunes the system through the zero-temperature critical point and might 
be related to the pressure or to chemical composition. The coupling strength uq is assumed 
to be constant near the critical point. The relation of 6q to the above mentioned control 
parameter r will be given below. The wavevector k is measured from the ordering wave 
vector Q. C,q is a microscopic length, presumably of the order of kp^. A microscopic energy 
scale is given by T which one expects to be of the order of Ep. The summation index n labels 
the Matsubara frequencies !)„ = 27m ksT. Momentum integrals are cut off at 7-, energy 
integrals obtained by analytically continuing Matsubara sums are cut off at rT(^o^)^~^- In 
a fermi liquid one expects the dimensionless cutoff parameters A and F to be of order unity. 
The dots indicate terms which are irrelevant at the critical point. 

In Ref. [^ the critical behavior of the model was shown to be determined by perturbative 
(in Uq) scaling about the Uq = (Gaussian) fixed point. The scaling procedure used in 
maps Z onto a new Z' identical in form to Z but with rescaled parameters T{b), S{b), u{b). 
The solution of the scaling equations can be written 

t{b) = t¥ (2.5) 

u{b) = m6^-('^+^) + 0{u^) (2.6) 

5{b) = b^ ho + 2u{N + 2) r^' dx e[2-('^+-)l- /(2)(te^-, 5(e"))| + Oiu") (2.7) 

tS^''\b) = ^N r^' dx e-('^+")"/(°)(te^",<^(e")) + Sf^ (2.8) 

Here t = ^ with T* = ^ and u = ^ with u* = ^. The analytical expressions for f^'^\x, y) 
and f^^\x,y) have been obtained in and read: 

f^°\x,y) = / — coth(^- )arctan( — — — ) 



{2-11^ Jo IT ^ 2x y + A'^' 

.±T -^k'-^coth{T^)aTctan{—-—) (2.9) 

71 Jo {27iY ^ 2x ' y + k^' ^ ' 



f'K^.y) = ^4^^:^ / - coth(^^^^ '- 



{2nY Jo 7r ^ 2x 'e'^ + {y + A^) 
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+ -r^ / 7 — ^ coth r — — - 2.10 

71 Jo {27iY ^ 2x 'T^ + {y + k^y ^ ' 

Note we have included the y-dependence of /*^^^; this was neglected in Ref. [Q. Its inclusion 
does not affect physical results in an important way, but makes the mathematical structure 
more clear. 

We turn now to the solution of these equations. For this purpose, we rewrite ( ^.71 ) 

5{h) = b\r{t, 6,u)-u A(t, r{t, 6, u), b)) + 0{u^) (2.11) 

with 



r{t,S,u)=S + 2u{N + 2) / dxe^^-^''+'^^'' [f^^\te'^,r{t,5,u) e^"") - f^^\0,0)] (2.12) 

Jo 

6 = 6o + 2{N + 2)uf^\0, 0)/[d + z-2] (2.13) 



and 



A{t,r{t,5,u),b) = 2{N + 2) / dxe^^-^'^+'^^'' f^^\te''',r{t,S,u)e^^) (2.14) 

Jlogfe 



We can neglect the 0{u^)-teTms in ( |2.11| ), because they do not affect the leading nonanalytic 



behavior in u. The expression ( ^^.12] ) is calculated in Appendix ^. The important result for 
our subsequent considerations in this Section is that we can write r(t, 6, u) as 

r(t,5,u) = 5 + {N + 2)ut'^^ <Pd^,(^-^^ (2.15) 



with 0rf^2 given by Eq. 



To calculate the free energy we use the scaling equation [13 



F{t,5,u) = b-^'^+'^ F{tb\5{b),u{b)) + tS^^\b) (2.16) 



with F{t,5,u) being the free energy of a system described by the model ( p.2| ), calculated 
perturbatively in the parameter u and measured in units of ksT*. Investigation of the scale 
(b-) dependent terms up to 0{u) in ( [^.l(j| ) shows that we can write ( |2.1(j| ) as a trajectory 



integral in phase space |l^ : 
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1/ /'(JXJ 

F{t,5,u) = j^N [ da;e-(^+^)^/(°)(te"^5(e^)) + C(M^) (2.17) 



V_ 

T^' Jo 

where the (9(M^)-terms do not contain the leading nonanalytic behavior in u. Inserting 
expression (|2.11| ) for 6{e^) and performing a variable change we get 






d+z-2 

ut '^ ^ f , r .. 

1 A l,-,e^ 

r V ti 



e^^' 



f2.18) 



It is therefore apparent that we can write the free energy (apart from non-singular contri- 
butions) as 

d+z-2 

F{t,6,u)=t'^f{^,^^^^) (2.19) 

f- r 

Due to the relation (|2.15|) and the definition of A{t,r{t,6,u),b), this scaling behavior is 
preserved when taking derivatives with respect to temperature, and we find the specific heat 
coefficient to obey (1^). It is shown in Appendix ^ that the dependence on the variable 



d+z-2 

— — - — leads to negligible corrections to the leading scaling behavior. The latter is given by 






and is evaluated in detail in Appendix |B[ 

In order to characterize the specific heat near the transition line tc{S) we calculate the 
discontinuity (jump) of the specific heat as it is predicted by mean field theory, using the 
scahng equation (|2.16|) for the free energy. Very close to the true transition, the model ( p.2|) 
can be scaled to a classical Gaussian model given by the Landau-Ginzburg functional 

exp(-»M,./p,,,,p,_5|,l, (,21) 

S{$] ^jd'x {i|W? + ^(6)1^1' + v(bm\'f] (2.22) 

with v{h) = u{b)t{b). Applying the usual mean field approximation yields the discontinuity 
of the specific heat. Analytical details of this calculation are given in Appendix ^ The 
result is Eq. ^78| . The correction to the specific heat jump due to Gaussian fluctuations is 



also evaluated in this Appendix. The result is Eq. (p.9|). Non-gaussian critical fluctuations 



become important when the Ginzburg criterion is violated ||T3|. The Ginzburg criterion may 



be obtained by demanding that the Gaussian fluctuation contribution to the speciflc heat 

d — A 

is small compared to the mean fleld speciflc heat jump. The result is utcV^^ <^ 1. This 
is violated for t — tc < tc "" in 3 dimensions. The expressions for the Ginzburg criterion 
and the width about the transition line where it is violated correct the results obtained in 
1^, which were in error by factors of t].'^ . As the critical temperature is small, mean fleld 
behavior characterizes the transition correctly except from a negligibly small region about 

tc. 



III. RESULTS FOR 3D 

We give the results for the speciflc heat coefficient (normalized per mol) and for the 
transition temperature in the physically relevant case of 3 dimensions in this section. 

The formula describing the speciflc heat in the whole Gaussian region of the phase 
diagram is Eq. ( |B3D [which can be simplifled to ( [B4{ ) for z = 2]. However, in order to 
compare experimental data to the prediction of this formula it is necessary to extract the 
system-dependent scale T* and the scale factor which relates r{T = 0) to experimentally 
accessible control parameters like pressure or concentration flrst. For this purpose it is useful 
to have expressions for 7 in the quantum (T <^ T*r) and classical (T ^ T*r) limits in the 
phase diagram. Fits of the data to the appropriate functional forms provide the scales which 
can then be inserted into (P3|), which must in general be evaluated numerically. 

The scale factors r* (which we have set equal to 1) and T* are of course arbitrary. Only 
dimensionless combinations of the coefficients we present are universal. However, there is a 
natural definition of the scale factor which gives an estimate of the important microscopic 
scales of the system. We have defined T* so that for d = 3 and z = 2 the leading low-T 
non-singular contribution to the Gaussian model specific heat is ^ji-§;^- For a density of 
order 1, ^q"^ r^ kp r^ (37r^)3 so the specific heat coefficient would be of the order of the 



leading fermi- liquid term in ( p. . 1|) . 

Furthermore, it is important to note that the specific heat calculated in this communi- 
cation is only the contribution arising from the existance of a quantum critical point. A 
metallic system has an additional background specific heat given by ( |1 . 1| ) . Depending on 
the dimensionality and the dynamical exponent, the leading behavior or corrections to the 
leading behavior in (|1.1| ) are changed due to critical fluctuations. 

Antiferromagnetism {z = 2): The leading term of 7 is still the normal fermi liquid 
contribution, a nonuniversal constant 70. Only corrections to this leading behavior are 
affected by the critical point. In the quantum regime (T <^ T*r) we find: 

^ ~^° T* 6 (r*)3 60 rl ^ ' 

Thus the coefficient of the first ((9(T^)) correction to the fermi liquid behavior diverges as 
r — > in the quantum regime. The classical limit (T ^ T*r) yields 

d<.„.o.,...2 keN^ 15iV I2\i „5,„. keN^ 3N (2\i „3, r 



with (^(|) = 1.34149 and C(|) = 2.61238. The term linear in r dominates corresponding 
terms in the non-critical part of the specific heat because it is enhanced by a factor T^2 at 
low temperatures. 

Clean Ferromagnet {z = 3): In the quantum regime both the leading (constant) 
term and the T^ log T-correction of the normal fermi liquid expression (|1 . 1| ) are enhanced 
by critical fluctuations: 

7— - - - - ^ £ •»« '- ^ ^ ? '°<^) - ■ ■ (3.3) 



^.„„.,,„3 ^ _^ A^ ,„^ (T^ _ Ma miM) i, + . . . (3.4) 



The classical regime shows non-fermi-liquid behavior: 

■bNa at ^^ /T\ _ keN. 
T* 97r ^^\f) (T*)^ 127r2v^ T 

with r(|)C(|) = 4.28742. 

Dirty Ferromagnet {z = 4): In this case, the leading behavior in the quantum regime 

is of fermi liquid type and enhanced by critical fluctuations, but the first correction is not 



ferini liquid like, and the coefficient of the latter diverges when approaching the critical 
point. 



^ T* 6 (T*)i 64 U/ 4^ ri 



(3.5) 



In the classical regime, the leading behavior is determined by critical fluctuations and is not 
of fermi liquid type: 



ii\/-n\ 
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dassicaLz=A 



kBNA N r(^)c( 



(T*)! 87r2 sin(|7r) 



,;.y,;j.-l_kBNA N r(|)C(^ 



(T*)3 327r2 sin(|7r) t! 



+ ... 



(3.6) 



with r(^)C(|) = 3.15612 and r(|)C(|) = 5.20628. 



Transition region: In order to characterize the region near the transition line Tc{5) we 
give the result for this curve in 3D 

-5 



rp rp* 



V^^,{N + 2)u_ 
the expression for the specific heat jump predicted by mean field theory 

kBNA 



(3.7) 



AC 



-N{N + 2yvl,uT. 



(T*)^(^+^)' 
and the correction to the mean field result due to Gaussian fluctuations 



(3.8) 



C^^'^^^'^ion ^ j^^j^^ 



N 



167r V^+1 



V3,,{N + 2)u 



T 



3 

2 


\T - Tel 


[ T, J 



(3.9) 



The numerical values for V^^ are 1)3 2 = 0.52109, ©3^3 = 0.39396, and 1)34^ = 0.35875. 



IV. COMPARISON TO EXPERIMENTAL DATA 



The discovery of superconductivity in so-called heavy-fermion systems |1[ resulted in 
extensive studies of these materials. Some of them undergo magnetic phase transitions at 
low temperatures, others display no long-range order. It has become evident, however, that 
several of the latter are near magnetic critical points as well. As an example, modest doping 
of CeCue |TB|,0 and UPt3 ||IB|,P^ leads to antiferromagnetism. The influence of pressure 



variation on the critical behavior has been studied as well [17,20|. The reason for the 
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proximity of the heavy-ferinion metals to a quantum critical point is the extreme sensitivity 
of the two competing energy scales (fermi [Kondo] temperature for the fermi liquid state and 
RKKY exchange energy for the magnetically ordered state) in these materials to changes 
in parameters such as the lattice constant. Therefore the heavy fermion compounds may 
be an appropriate experimental system to study the implications of the theory presented in 
this communication. 

Application to CeCu2Si2 

Measurements of the temperature-dependence of the specific heat have been carried out 
at compounds Cei-^M^CuaSiz with M=La,Y,Lu and x = 0.0,0.05,0.1,0.2,0.5 and 0.9 [g. 
Doping with La expands the lattice, whereas doping with either Y or Lu contracts it. As 
the data for La-doped samples show a rather complicated behavior, we will not comment 
on them in this work. Diluting CeCu2Si2 with Y,Lu indicates, however, that this material 
is near a quantum-critical point. Based on the assumption that this is an antiferromagnetic 
instability, we present an interpretation of the data for Cei_a;Lu2:Cu2Si2. The numbers and 
conclusions which can be drawn from the experiments carried out with Y-doped samples do 
not differ significantly from those given below. 

In Fig.|^ the data for various concentrations x are shown. The parameter tuning the 
system through the quantum critical point 6 = r{T = 0) is defined as 6 = ^^^. If the 
pure compound is right at the critical point [xc = 0), this system would be in the classical 
regime and we would expect to find the low-temperature specific heat behaving according 
to ( |3.2|) . Indeed the low-temperature (T < 5K) data for x = fit nicely to the square 
root temperature-dependence. This fit yields the temperature scale T* = 5K which is of 



the order of the Kondo temperature in this material: T^ = lOK |^1|. Furthermore, the 
first correction to the leading behavior in the classical regime is linear in r and therefore 
linear in 6 as well. As we will check below, the compounds for x = 0.05, 0.1 are still in the 
classical regime. We take the measured values for ^(T = I.IK) given in [Q at concentrations 
X = 0.0 and x = 0.05 in order to determine the scale x*. We find x* = 0.3, hence for the 
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quantum-classical crossover temperature T^ross = ^ ^ = 17x which is reached right at 
the beginning of the measurements for x = 0.05,0.1, whereas for x = 0.2 the quantum- 
classical crossover is accessed experimentally. The data for x = 0.5, 0.9 and T > 5K cannot 
be explained within the framework of the theory presented here because the systems are 
already too far away from the critical point. Taking the values for the nonuniversal scale 
factors T* = 5K, x* = 0.3, we can apply ( |B4| ) to predict the temperature dependence of the 
specific heat for x = 0.05, 0.1, 0.2. As Fig.|^ shows, theory and data agree quite well for low 
temperatures. Please note that we neglected any effect of the dangerously irrelevant variable 
u. The coincidence for x = 0.2 at low temperatures is particularly nice because there the 
system is in the quantum regime, whereas we extracted all information for the theory from 
the classical regime. There is apparently a systematic deviation from the theory occurring 
at lower temperatures for larger x. We do not understand this deviation at present. We 
suspect that the deviations (starting at lower temperatures for increasing concentration) are 
due to noncritical fluctuations which become important farther away from the critical point, 
but the observed interplay between T and x is not expected in our model. 

V. CONCLUSION 

We have determined the universal functions needed to quantitatively analyze the behavior 
of the specific heat near a quantum critical point with dynamical exponent z = 2,3,4 in 
d = 3 spatial dimensions and illustrated, by an example, the use of the equations. The 
results are summarized in Eqs. ( p.l|) -( ^^ , (|B^ ) and (p^, respectively. Depending on the 
dynamical exponent z, the leading or subleading behavior of the normal fermi liquid specific 
heat coefficient is affected due to the proximity to the quantum critical point. 

The physical properties of these transitions are governed by four non-universal scales: a 
temperature scale T*, a scale r* for the Hamiltonian parameter tuning the T = transition, 
a dangerously irrelevant variable u, and a length scale ^o- Once these are known, all the 
behavior is fixed. The specific heat does not involve the length scale ^o but provides sufficient 
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information to overdetermine the remaining parameters. 

An important open theoretical problem is to understand the nature and magnitude of 
the corrections to the universal behavior. On the experimental side, it would be useful to 
study in detail the specific heat of a system to which the theory applies. 
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APPENDIX A: QUANTUM PARAMETER AND CRITICAL TEMPERATURE 



In this appendix we will calculate (|2.12| ) and obtain an expression for the critical tem- 
perature. 

We first extract the explicit temperature dependence of r when writing 



r(t, 6, u) = r{r(t, 6, u),u) + r(t, r{t, 6, u),u) 
f{6, u) = 6 + 2u{N + 2) dx eP"('^+-)l- /(2) (o, r (t, 5, u) e^") - Z^^) (0, 0) 



(Al) 

(A2) 

f{t,6,u) = 2u{N + 2) /" da; e !'-(''+* f/(2)(te"",r(t, 5,m) e'") - /(')(0, r(t, 5, u) e"")] (A3) 

JO L J 

Further inspection of Eq. (|A2| ) reveals that 

/ 

f{5,u) = 5 l + 2n(A^ + 2)- 



d,f'\^^v)U=,\ . . ... . .. I Oir'-^) ifrf + z-6^0 



d + z-A 



2u{N + 2) 



0{r^\ogr) iid + z-Q = Q 



(A4) 



where we omitted 0{u^)-ieiins, like those we have already neglected when writing ( p.ll ). 
The dangerously irrelevant variable u leads to a rescaling of 5 by a u-dependent factor, 
which can be absorbed into the scale (unit of measurement) of 5. Note that the nonanalytic 
contributions in r due to the variable u are nonleading and therefore negligible near the 
critical point. 

Neglecting nonleading (analytic) contributions we can write for (A^) 



d+z-2 



with 



fit, S, u) = u{N + 2) t^^ <Pa,z ( ^ ) + Oiut^' 



Kd ^ [^ . ^ [^ . e(cothe-l) 
W(a^) = /^ w+i 8 / dyy ^ / de ■ 



(A5) 



(A6) 



(27r)^+i Jo Jq Ae'^y^ + {xy + lf 

The two limits -^ -C 1 and ^ ^ 1 correspond to the classical and quantum regimes, 
respectively 0. We find 



^^,2(2;) ^ Vd^z for a; — > (classical) 
V^(^) ^ ^ ^ for a; —* 00 (quantum) 
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with 






i^-nY ^ sin 



-TT 



22 

Numerical values for P^ ^ in three dimensions and for the dynamical exponents considered in 
this communication are T'3_2 = 0.52109, "Ds^s = 0.39396, and P3 4 = 0.35875. The term with 
explicit temperature dependence in the quantum regime is t^r 2 which is small compared 
to the contribution r 2 to f which we have already neglected. The general result is 
therefore 

r (t, 5, n) = 5 + (iV + 2) M r^^ 0d,, (^^) (A8) 

= (5 in the quantum regime (A9) 

= (5 + P(i^^ (A^ + 2) Mt ^ in the classical regime (AlO) 

The true transition occures when r vanishes. Therefore we find an expression for the 
transition temperature tc by setting ( |A10| ) equal to zero and solving for t: 

-6 



tc 



Pd,. {N + 2) u_ 



d+z-2 

(All) 



APPENDIX B: QUANTUM-CLASSICAL CROSSOVER 



In this section we investigate the expression (|2.20|) for the leading contributions to the 
specific heat coefficient, which arise from the quantum-classical crossover. 
The expression ( p.2(]| ) can be written 



^'''"''''' e,{2T^Y+'J^ Jo sinhV(2te)2a;^ + (rx + l)2 ^""'^ 

-" eo{2nr^'Jo """sinh^J^ {x + 1)^ + ^ x^ ^ ' 

In the case oi d ^ z, the dependence on the cut-off A gives rise to analytic terms 0{r,t^) 
only, which we will neglect in the following. To discuss the case d = z, however, the existance 
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of a lower boundary of the x-integral is essential for technical reasons because logarithms 
occure. But even in this case, the nonanalytic behavior of the heat capacity does not depend 
on the cut-off. 

For (i = 3, z = 2 it is possible to perform the x-integrals and we find for the singular part 
of 7 which arises from the quantum-classical crossover: 



In the next paragraph we study the asymptotic behavior of jieading in the classical and 
quantum regimes corresponding to the limits ^ <C 1 and ^ ^ 1, respectively. For sake of 
simplicity, we do not write the factor t? {•^'y+i explicitly. Results of the following paragraphs 
have to be multiplied by this term to obtain the correct results. 

Classical regime: To investigate the classical limit it is most convenient to use Eq. 
( p3D and expand in the parameter -^ (which is small in the classical regime). The leading 



behavior is 



classical4l^z _ 27r r(^ +2)C(^ + 1) it£ ,^^. 

1 leading " ~ '. /^.^ \ ^ + ■ ■ ■ l,DOj 



.in (^vr) 
in the case d ^ z^ whereas for d = z we have 



47r2 
classicaLd=z ^" i , . /nn\ 

lleading " ~^ °^ ^ ' ' ' ^ ^' 

The first correction to the leading term is obtained by an expansion in the small parameter 
and reads 

vr r(^ + 2)C(^ + l) r , , 



2z sin (^^vr) t 



2 + z-d 



Quantum regime: In this limit we use an expansion of (|B2|) in the small parameter 
^. The leading term for d ^ z is 

iZZT'^' = r" 'Sniz - rf) r^ + . . . (B8) 

In the case d = z we find 
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quantum, d=z 
Heading 



'r i°sa^+--- 



(B9) 



For z = 2 the corrections to the leading behavior are obtained by expanding (B^) in —^ 



This yields the asymptotic series 



quantum, djtz=2 
Heading 



-47rV^ ^(-l)"(27r)"|5; 



2n+2| 



n=0 



r(2n - \) n_ 
r(-i)r(2n + i) \r 



2n 



2n^ 



1 

f 2 



15 ra 



(BIO) 
(Bll) 



where Bm denote Bernoulli numbers. Such an expansion is inappropriate to find corrections 
in the cases z = 3,4 because these are nonanalytic in the quantity ^. It is more convenient to 
rewrite ( p3D for large ^ and subtract the already calculated leading term. The corrections 
are then given by 



d I -I d 



z — 2 ^z-2 v/o 



de 



5Z-2 



o+l 



dx 



X 2(^-2) 



sinh^ e Jf2i4ll 1 + X 



which gives for d = z = 3 



(B12) 



and d = 3,z = 4 



15" ;:^i°g^ + --- 



15 /7r\l ^,5, t^ 



(B13) 



(B14) 



APPENDIX C: DANGEROUSLY IRRELEVANT VARIABLE AND THE 
QUANTUM-CLASSICAL CROSSOVER 



In this appendix we investigate the dependence of (|2.18|) on the argument 



d+z-2 

"* " in the 



quantum and classical regimes. We show that it leads to corrections to the leading behavior 
given by (|1.2|), which is calculated in Appendix 0. 

In the quantum regime (t -C r^/^) we can write (|2.18| ) as 



^0 "'2iogr vr ' 



d+z-4 

1 — ur 2 A 



1 e^ 

,2/2' ''^ 



,2x 



(CI) 
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As the parameter -^ is small, an expansion in this quantity is reasonable. The term ur 2 
is small compared to unity near the critical point {u is assumed to be of order unity, and r is 
small), it leads to corrections to coefficients in the leading behavior in the quantum regime, 
the latter is given by an expansion of (P^) in the small parameter -^. 



Similarly, in the classical regime we can expand (|2.18|) in the small quantity -^. The 



d+z-2 



d+z-4 



correction due to the argument — — - — is here ut"^"^ " , which again has to be compared 
to unity. Near the critical point, this correction is therefore negligible. However, a careful 
calculation of the terms arising when taking derivatives w.r.t. temperature shows that a 
nonanalyticity arises due to the temperature dependence of r{t, 6, u) in the classical regime. 
This leads to a divergence of the specific heat which is calculated in Appendix |^. 

In conclusion we find the variable u to give rise to corrections to coefficents in the limiting 
expansions of 'pleading (cf. ( p2D and (p3|), respectively) only - apart from a narrow region 
about the transition line tc{5), where the Gaussian correction to the discontinuity of the 
specific heat becomes important (see discussion in the last paragraph of Section ^). The 
coefficients in the expansion of 'pleading are of order unity, whereas the corrections in ur 2 
and ut s are small near the critical point. 

APPENDIX D: SPECIFIC HEAT JUMP AND DIVERGENCE NEAR THE TRUE 

TRANSITION 



The characteristic behavior of the model ( p.2|) near the transition is calculated in this 
section. As it has been shown in |Q, (|2.2|) scales to a classical theory given by ( |2.22|) . For 



this theory, the discontinuity of the free energy reads [13 



Using the relation -F(l) = 6^*^'^+^-' F{h) in order to match the scaled model to the original 
one, we find the expression for the discontinuity of the free energy of model ( |2.2| ). We insert 
( [^.6D for u{h) and substitute 5{h) = b"^ r(t, 6, u) with r(t, 5, u) given in Eq. ( JAIOD . The latter is 



justified according to ( p.ll| ), noting the fact tliat tlie quantity A(t,S,u) becomes negligible 



in the classical regime. Differentiating twice w.r.t. temperature and taking into account 
that r{t,6,u) vanishes at the transition, one ends up with the quantitative expression for 
the specific heat jump: 

AC = Vl^N{N + 2yutl'^^^''~^'^ (D2) 



with the numerical constant ^^,2 given by Eq. (|A^). Please note the anomalous exponent 
of tc in this expression. In theories where the distance to the critical point is measured by 
{t — tc) - classical Landau theories - we find that AC is proportional to tc- For our case of 
a quantum critical point, the specific heat jump turns out to be smaller because additional 
(in general noninteger) powers of tc occure. 

Although the effective dimensionality of the theory ( p.2|) is d + z (which is greater than 4 
in the cases discussed in this communication), a divergence of the specific heat occures very 
near the true transition (for vanishing r{t,6,u)). This is due to the fact that the quantum 
fluctuations become irrelevant and a classical description of the transition becomes more and 
more appropriate when the critical temperature is approached. To calculate the divergence. 



which has its origin in the temperature dependence of r{t,6,u), we start from ( p.l8| ) and 
express the temperature derivatives (which we have to take in order to obtain 7) in terms 
of derivatives w.r.t. r(t, 6, u). Apart from contributions which are negligible corrections like 
those calculated in Appendix |C|, and analytic terms, we get 

transition . ^ .2 ^^^N ^,i±^ f-- , ^^4^^^, foo £ 8 {y + 1) 



^trans^t^on ^ ^Q^^y a ^^ , ^^^_ 1 , ^^ ^^^^(^ j) 



(27r)'^+i h '' Jo -2^y^'[e' + {y + irf 

(D3) 

The divergence occures for r{t, 6, u) -^ 0, we can therefore expand the hyperbolic cotangent 
in its argument. Only the leading term gives rise to a divergence in r^transition ^ which is 

7*™*^°" = (5tr)^i^tr^ rdyi-^-'iy + l) H de ^ ^ , /, ^,,,, (D4) 



(27r)'^+i Jo ^ ^ ''^ ' Jo [e^ + {y + l)^]^ 

Performing the (elementary) integrals and using (|A10| ) to evaluate dtr yields finally: 
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c 



transition 



NKd{d-2) fd + z-2' 



d-A 

2 



(27r)'^-il6sin^7r V z 
The numerical constant P^^ ^ is given by Eq. (| 



Vd,,{N + 2)utc 



d+z-2 



-t-L. 



(D5) 
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FIGURES 
FIG. 1. Phase diagram for a quantum-critical point which is controled by a zero temperature 

parameter 6. The shaded region is the (magnetically) ordered phase. The solid line denotes the 

transition curve Tc{6), the dotted line corresponds to the quantum-classical crossover which occures 

when T becomes of the order of r^'^. 

FIG. 2. Comparison of experimental data measured from CeCu2Si2 to the theoretical prediction 
for d = 3,z = 2. The temperature scale T* = 5K has been extracted from the data for x = 0, 
assuming the pure compound being at the critical point. The scale x* = 0.3 is determined by the 
knowledge of T* and the datapoint at z = 0.05, T = I.IK. Solid curves represent the theoretical 
prediction, symbols are data points. 
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